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The generalized Langevin equation revisited: Analytical expressions for the 
persistence dynamics of a viscous fluid under a time dependent external force. 


Wilmer Olivares-Rivas^, Pedro J. Colmenares^’* 

Grupo de Quimica Tedrica: Quiraicofisica de Fluidos y Fendmenos Interfaciales (QUIFFIS). Departamento de Qmmica - Universidad de 

Los Andes. Merida 5101, Venezuela 


Abstract 

The non-static generalized Langevin equation and its corresponding Fokker-Planck equation for the position of a viscous 
fluid particle were solved in closed form for a time dependent external force. Its solution for a constant external force was 
obtained analytically. The non-Markovian stochastic differential equation, associated to the dynamics of the position 
under a colored noise, was then applied to the description of the dynamics and persistence time of particles constrained 
within absorbing barriers. Comparisons with molecular dynamics were very satisfactory. 
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1. Introduction. 

The dynamics of particles in a potential field is usually 
investigated by molecular dynamics (MD). However, since 
the time scale of such a procedure is of the order of fem¬ 
toseconds, the algorithms are in general very expensive in 
computer consuming time. To overcome this difficulty, one 
generally appeals to mesoscopic descriptions on which the 
time scale is larger. It gives a glimpse about the dynamical 
behavior of the system at the time scale of the technique. 
There are three levels of description. The most commonly 
used description is the so called Langevin equation (LE), 
where the dynamics of a Brownian particle in phase space 
is described by the Markovian set of stochastic differential 
equations (SDE)[iHi 

= -l{z)v{z,t) + r]{z)£,{t)+F{z,t), 

Z = V, (1) 

where z{t) and v{z,t) are the position and velocity of the 
tagged particle of mass m at time t, respectively; 7 ( 2 ) 
is the position-dependent phenomenological fluid friction 
coefficient; F{z,t) is the external force field, r]{z) is the 
intensity of the stochastic force and ^{t) is a zero-mean 
Gaussian white noise term. The static LE can be solved by 
a double integration or alternatively, throi^h the solution 
of a second order differential equation[l[ Q. In fact, an 
analytical solution for a confined fluid has been previously 
given 0 . 
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An important simplification of the LE with a wide range 
of practical applications is the so-called Smoluchowski 
limit also known as the high friction limit (HFL). The 
most simple derivation consists on assuming that the fluid 
has a large enough friction coefficient 7 so that the velocity 
in the LE shows a fast relaxation to a quasi equilibrium 
[^, with dv/dt = 0. Thus, position evolves with time ac¬ 
cording to 

-=vAz) + V^om, ( 2 ) 

where Vp{z) = F{z,t)/j{z) is the induced velocity drift 
and Do = {l/2){ri{z)/j{z))‘^ is interpreted as a diffu¬ 
sion coefficient. The corresponding master equation for 
the probability density corresponding to this SDE is the 
Smoluchowski equation (SE). Clearly this is only valid 
in the long time limit, also known as the Smoluchowski 
regime. The third stochastic description is an improve¬ 
ment without restriction on the time scale, known as the 
non-static or generalized Langevin equation (GLE). It is 
obtained modifying the dissipation or friction term 

^ dv{z, t) _ f Y(^z,t — t')v{z,t') dt' 

at Jo 

+Rit) + Fiz,t), (3) 

where the kernel r(z, t — t') expresses the memory or retar¬ 
dation effect on the movement of the fluid particle due to 
the collective hydrodynamic response of the surrounding 
fluid 0 and, R{t) = r]{z,t)J(t) a colored-Gaussian fluctu¬ 
ating driving force. Now the random noise function J (t) is 
not a white noise. If the retardation is omitted, the kernel 
is static r(z, t—t') = a{z)6{t—t') with. a{z) = 7 ( 2 )/m,and 
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the LE, Eq. o, is recovered. This is equivalent to take the 
noise as being white, namely, J(t) = ^{t). However, when 
inertial effects are not negligible -the friction coefficient is 
not dominant in the dynamics- the position and velocity 
are driven by non-markovian random terms having a finite 
correlation time. This is the situation for real non homoge¬ 
neous fluids and systems with strong boundary conditions 
or in the presence of chemical reactions. 

A well known work obtaining a master equation for the 
probability density | Co) for a general C process, as 

a Kramers-Moyal[7j cumulant expansion, is due to Hanggi 
in 1978[ii]. However, the GLE, Eq. (Ip . was discussed 
for a viscous fluid by Chow and Hermans [ifl as far back as 
1972 and a master equation for that generalized version, 
was nicely derived and solved by Dufty[Il| in 1974, being 
commonly referred to as the Fokker-Planck equation. In 
the absence of external forces it was extensively studied 
by Adelmanfl^. Fox[l^|. Volkov and Pokrovsky [l^ and 
later by Rodriguez-Salinas [l5j |. Numerical approaches for 
general first-order SDEs with colored random forces have 
been also presented[l^ [I^. Budini and Caceres[l^, nu¬ 
merically obtained the velocity distribution associated to 
the GLE SDE for arbitrary noise and memory kernels, but 
with no external force. They found that the interplay of 
noise structure and dissipation is an important issue to 
consider, in order to achieve the stationary steady state of 
the probability density. 

In this paper we analyze the dynamical or general¬ 
ized Langevin approach, Eq.(|31), for the motion of an 
interacting fluid particle in an external time-dependent 
field. We obtain the corresponding generalized Fokker- 
Planck equations (GFPE) for the Bayesian probabilities 
p{v,t\vo,0) and p(z, t|zn , un, 0). Then, following Chan¬ 
drasekhar’s approach 1 18 j . we obtain the vq averaged prob¬ 
ability densities p{v,t) and p(z,t\zo,0) giving analytical 
formulas for their moments. We shall discuss the fact that 
this equation can be written as a much simpler statisti¬ 
cally equivalent Markovian first-order SDE. Revisiting the 
master equation associated to GLE, we get consistent an¬ 
alytical results for the survival probability in a viscous 
fluid and evaluate important dynamical properties as the 
mean square displacement (MSD) for a fluid bounded by 
absorbing barriers. 

Analytical expressions for Chandrasekhar conditional 
probability density p(z,t\zo,vo) and the corresponding 
Rayleigh type GFPE in the presence of an external force 
are presented in Section[21 That is, with the method based 
on the solution of the stochastic Liouville eauation (l9l - [^ . 
we find a well defined GRFPE, whose analytical solution 
matches Chandrasekhar’s lemma [l^. 

Section [3] deals with the GLE’s velocity and position 
space distribution moments for any given external time 
dependent force F{t). First, we revise the GLE z-space 
to find the moments of the probability density p{z^t\zo) 
averaged over the initial velocity, unifying all the equiv¬ 
alent master equations and SDE. The known results on 
the velocity space probability BP are complemented 


to include a time dependent external force. 

Considering an exponential decaying friction coefficient 
kernel and an external constant force, we derive in section 
|4j analytical results for the moments of the distributions 
in z and v spaces. The application of the GLE theory to 
the problem of survival probability and first passage time 
of a constrained fluid is presented in section [S] 

Finally, we include two appendices to show an unified 
view of the velocity-fluctuation coefficient in the GLE ap¬ 
proach and the classical LE and SE results. 


2. The position generalized Rayleigh-Fokker- 
Planck equation (z-GRFPE). 


In this section, we use the simple procedures well dis¬ 
cussed by AdelmanB- FoxB Sancho et al. 211, to 
obtain the GFPE type master equations associated to the 
position of the particles of a viscous fluid satisfying the 
GLE. It has been shown [3, Bl fFat even for strongly non 
homogeneous fluid the application of molecular dynamics 
and LE equation can be carried out by sampling the z 
space in virtual layers where the local mean force is as¬ 
sumed piece-wise constant. Thus, we shall assume that 
the time dependent force is a constant as a function of po¬ 
sition, namely F{z,t) = F{t). Besides the fact that even 
the LE has no solution for a general F{z, t), it will greatly 
reduce the mathematical complexity of the equations. 

Following the standard method used by Fox[Il| , we ap¬ 
ply the k-Laplace transformation to both sides of the GLE 
differential equation, Eq. and find after collecting 

terms, an expression for the Laplace transform of the ve¬ 
locity v{k) 


= Xv{k) 


v„ + 


■[m+Fik)) 


(4) 


where is the initial velocity. The function Xv{k) is 


Xvik) 


1 

k+f(ky 


(5) 


with F(fc), R{k)^ and F'(fc), being the Laplace transforms 
of the memory kernel, the colored-noise internal force, and 
the external force, respectively. 

The fundamental Green function Xv{t) is the inverse 
transform of Xv{k). From the definition in Eq. @ we 
have the useful Volterra relationship 


Xv{t) = = - j Xv{s)T{t-s)ds. (6) 

Before solving Eq. o, it is useful to define two related 
functions. So, the integral of Xv{i) is denoted as Xz(t) and 
in turn, the integral of Xzit) is denoted simply as x{k), 
without any subscript 


Xzit) = / Xv{s)ds, 

Jo 


(7) 


2 







xit) = 

/ Xz[s)ds 
^0 

(8) 

= 

/ {t- s)xv{s)ds. 

^0 

(9) 

Clearly Xz = 

Xv and X = Xzi with initial 

conditions 


Xv(0) = 1, Xz(0) = 0, and x(0) = 0. 

Now, inverting both sides of Eq. we get the particle 
velocity v(t) = z(t) and by integration, its position z(t) 

v(t) = v{t) + ipy{t), (10) 

z{t) = + (11) 

The colored noise force response functions (fiy (t) and (pz (t) 
in the velocity and position components, are defined re- 


spectively as: 





1 

/■* 


Py{t) = 


/ Xv{t- s)R{s)ds, 

(12) 


m , 

Jo 



1 

/■* 


‘^z{t) = 

m , 

/ Xz(t- s)R{s)ds. 

Jo 

(13) 


The drift components of the velocity and position are 
actually the averages over the noise distribution, namely, 
v{t) = {v{t))^ and ^{1) = {z{t))^. They are given as 


vit) = VgXv{t) + (14) 

z{t) = zo + v„Xzit) + (l)z(t). (15) 


Here vq and zq are the initial velocity and position at t = 
0, and the corresponding velocity and position response 
functions to the external force, (^v(t) and ipzit), are given 
by 


(j)y{t) = 

1 

[ Xv{t- s)F{s)ds, 

(16) 

m , 

Jo 

4>z(t) = 

1 

m . 

[ Xz(t- s)F{s)ds 

Jo 

(17) 


4>vis) ds. 



From the solution for the velocity, Eqs. (fTUI) . (Upland (fTil) . 
it can be shown that the susceptibility Xv(t) is directly 
related to the velocity auto-correlation function (VAC), 
Cy{t) = {v{0)v{t))j^, namely, Xv{t) = Cy{t)/cy{0). The 
noise and external force response functions also have the 
simple relationships Pzit) = (pv(t) and (fzit) = (t>v{t). 
They are functionals of the susceptibility function Xv{i), 
namely (pv{t) = Pv{Xvit)]- Consequently, the conditional 
probability and the position itself are functionals of Xv{i), 
i.e., p(z,t|Zo,'Uo,0) =p[xv{t)\ and z{f) = z[xv{t)]. 

The mathematical and physical consistency of the prob¬ 
lem requires to know the statistical properties of pv (t) and 
(pz{t). Above all, they should have zero mean averages, 
that is {(pv{t))R = 0 and therefore {(pz{t))R = 0. They are 
described in appendix A, where in particular, we derive 
the appropriate fluctuation-dissipation theorem. 


The velocity space has been amply studied. So, in this 
work we shall focus on the less studied position space. Nev¬ 
ertheless, in Appendix A, besides summarizing the proper¬ 
ties of the color noise functions, we discuss the equivalent 
Fokker-Planck type generalized master equations for the 
velocity conditional probability p{v, f | uq, 0) associated to 
the GLE, Eq. ©. 

In order to derive the corresponding FPE for position, 
let p{z, t I zo, uo) be the conditional probability distribution 
of finding the particle, say at position z at time t, given it 
started to diffuse at z^ with velocity uq at t = 0. To obtain 
the FPE for the evolution of p(z,t | zo,uo), associated to 
Eq. @, we will use a method originally dev elop ed by 
Sancho et al. 211 and applied in many problems |23l l24l| . 

For a given realization of the noise ipv{t), Eq. (fTOl) de¬ 
scribes a flow in z-space. The density of this flow evolves 
in time according to the stochastic Liouville equation 


df{z[ipy],t) 

dt 


dz 


f{z[Pv\,t) 


dz[py 


dt 


(18) 


where f{z[(py],t) is the probability density of the flow. 
Taking into account all realizations o f iPy (t), the Liouville 
equation turns into an ordinary SDE[19|. As pointed out 
by van Kampen the probability density of the fixed 

realization z at time t, p(z,t | zo,uq), can be obtained by 
averaging the function f{z[(pv],t) over the distribution of 
the colored noise, namely 


p{z,t I zo,'yo) = {f{z[pv\,t) 


(19) 


Then, by replacing Eq. m into Eq. m, the equation 
satisfied hy p(z,t \ zo,vo) reads: 


dp{z,t I zo,uo) 
dt 


= -v{t) 

dz \ 


dp{z,t I zo,uo) 
dz 

(f{z[ipy\,t) (Py{t) 


( 20 ) 


Here, the subindex tpy is a remainder that the probability 
density of the colored noise has to be employed in the 
calculation of the average. Since ipy(t) is a zero-mean 
Gaussian noise, the cross correlation {f{z[ipy],t) (/?„(t))_, is 


given by the formula of differentiation due to Furutzu 261 
Novikov [23 and Donskerf^ 


- {f{z[^v],t) Py{t)'j = 


tp. 


f ^\ / 9f{z[pyft) 5z[py\\ ^ 

( 21 ) 


where the functional derivative dz^Py]/5py[s) = I was ob¬ 
tained from the definitions, Eqs. (El) and m, and van 
Kampen definition, Eq. (HU), was used. 

Then, carrying out the proper substitutions in Eq. (1201) . 
we finally find that the probability density satisfies the 
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following partial differential equation 


/dp{z,t\zo,vo) 

V dt 


dp{z,t\zo,vo) 

-’W—&— 

n d‘^p{z,t\zQ,VQ) 


( 22 ) 


where the time dependent diffusion term Dq(t) was written 
as 



(23) 


in which C{t, s) = {pv{t) fv{s))ji is the two times correla¬ 
tion of the colored noise response function. This differen¬ 
tial equation is often referred as a Fokker-Planck equation. 
It is actually a Kramers-Moyal z-space master equation, 
analogous to the Rayleigh equation in v-space. To empha¬ 
size the vq dependence and the fact that the associated 
SDEs, Eqs. (ITUl) and (fTTl) . constitute the exact solution 
for the GLE, we shall refer to Eq. (0^ . together with the 
z-moments derivatives, Eqs. (HI and (1231), as the z-space 
generalized Rayleigh-Fokker-Planck equation (z-GRFPE). 
We reserve the acronym GFPE for the master equation 
associated to the vq averaged p{z, t \ zq) to be obtained in 
next section. 


3. Time dependent external force: Alternative 
SDE GFPE for position and velocity. 

In this section we write equivalent forms of the z- 
GRFPE, Eq. (1^ . corresponding all to the non-static 
GLE, Eq. ([3]), and give the analytical solutions ioi p{z,t \ 
zo, Vo) and its initial velocity average p(z, t\ zq) = {p{z, t \ 
Zo,vo))voj for any given retardation kernel F(t) and exter¬ 
nal force F{t). 


3.1. Solution of Position GFPE 

First of all, applying the simple linear transformation 
q{t) = z{t) — z(<), with the variables J{t) and v{t) = z{t) 
as defined in previous section, Eqs. (TT^ a nd (TTKll . Eq. 
(EH), reduces to a diffusion-like equation [29| 

(dp{qG\qo,vo)\ ,_,.,d'^p{q,t\qo,vo) 

V dt dq^ ’ ^ ’ 


Its solution with an initial condition S{q — qq) is a Gaussian 
centered at q^ = 0. In terms of the original variables, the 
conditional probability density for a given realization z, 
starting from z^, given an initial velocity vq in a GLE 
process described by Eq. (El) is then 


p{z,t I zoWo) 



2 ^ 2 (f) . ’ 


(25) 


where the standard deviation aq{t) = {{z — z{t))"^)^ is de- 
hned as 


(Jqit)= 2 / Dq{a)da. 


(26) 


4 


The correlation C{t, a) is evaluated in Appendix A, Eq. 
(IA-241) . Using it into the expression obtained above for 
the function Dq{t), Eq. (1331) . we can readily evaluate 


D,{t) = 


Cit, s)ds, 


KT 


m L 


Xv{\ t 


s\)ds-Xv{t) Xvis)ds 

Jo 


4>v{t) / (/)«(s)ds, 

Jo 


(27) 


Using the definitions of the functions Xz (t) and (fz (<) 


D.it) = ^Xzit) 


i-x«(f) 


- 4'v{t)4>z{t). 


(28) 


A second integration gives the second moment of g = 
[z{t) — 'z{t)\ in terms of the function x(t) 


(Zq{t) = 


k^T 


2 x(f) - xl{t) 




(29) 


The probability p{z,t \ zq), of finding the particle at po¬ 
sition z{t) at time t given that it was initially located at 
Zo, independently of the initial velocity, is obtained aver¬ 
aging the conditional probability p{z,t \ zo,vo) over the v„ 
Maxwellian distribution 


Pi{vo) 


1 



exp 


mvo 

2^ 


(30) 


A simple integration of the Gaussian in Eq. (1251) gives 

/ OO 

p(z,t I zo,vo)pi{vo)dvo 

-OO 


\/27rcr^(t) 


exp 




2al(t) 


where 


(Zzit) = 

= 2 




q 

KT 

m 


m 


xW - 


(31) 

(32) 

(33) 


Therefore, the diffusion coefficient for p{z,t \ Zq) is: 

Dz{t) = Dq{t) -h ^^Xzit)Xvit) 

m 

= ^^^Xzif) - K{i)4>z{t). (34) 

m 

We should point out the fact that, since the particle is free 
to move in the entire space, p(z, t \ zg) is referred to as the 
unbounded position probability density in one dimension. 
It is in fact the solution of the exact generalized Fokker- 
Planck equation for position z (z-GFPE) 


' dp{z,t I Zq) 

^ dt 


) = 

/ Z 

+ T>z{t) 


dp{z,t |zo) 
dz 

d'^p{z,t |zo) 

az 2 


(35) 








































with the boundary and initial conditions 
p{-oo,t\zo) = p{QO,t \ Zq) = 0, 


piz,t = 

0 1 Zo) = 5{z-zq), 

(36) 

where 


II 

o 

II 

(37) 


= ZQ + (j)z{t), 

(38) 

zo 

= {zit = 0)U- 

(39) 


For an unbounded particle, using these last relation¬ 
ships, it is straightforward to write the MSD with respect 
to the initial position, denoted as cr^(t) without any sub¬ 
script, i.e., 

a^(t) = MSD{t) = (([z(t) - zo]")„).o, (40) 

we get 

+ <(>^(0 = 2—x(t). (41) 

m 

This result can also be obtained directly from Kubo’s vq- 
averaged VAC, c„(t) = {k^T/m)xv{t) 


we get the velocity-independent generalized diffusion 
equation (GDE) 

l' dp{Z,t\Zo) \ _ d^p{Z,t |Zo) 

V dt Jz(t) ^ dZ2 ’ ^ ^ 

with a time dependent diffusion term Dz{t) given by Eq. 
((Ml) . According to this, it is identical to D{t) only for 
zero external force. In that case, the common Sutherland- 
Einstein limit Dz{t) = D{t) —>■ Dq is obtained in the 
Smoluchowski limit where axz(t) —> 1. 

This generalized unbounded diffusion equation satisfies 
the initial and boundaries conditions: Zq = 0 and 

p{Z,t = 0\Zo = 0) = d{Z), 

p{-oo,t \ Zq) = p{oo,t \ Zq) = 0. (47) 

3.2. Solution of velocity GFPE 
The z-space master equation, Eq. (EH), gives the prob¬ 
ability density p{z, t\zo, vq). It is associated to the process 
defined by the GLE, with v{t) = z{t) given by Eq. (ITUl) . 
which can be written as 

v{t)=VQXv{t)+‘p{t), (48) 

with the total internal plus external force response func¬ 
tion defined as 


cr^(t) = 2 f (t — s)c„(s)ds. (42) 

Jo 

The regular definition of the diffusion coefficient commonly 
related to experimental data or molecular simulations is 

D{t) = = [ c^{s)ds = —Xz{t). (43) 

2 dt Jq m 

It is interesting that D{t) so defined is not explicitly depen¬ 
dent of the external force, while the coefficient directly as¬ 
sociated to p{z,t I Zq), namely Dz{t) = D{t) — (f)y{t)(j)z{t), 
does depend on F{t). It is also common to find the dif¬ 
fusion coefficient defined as an extension of the SE result 
(T^(t) = 2Dt, i.e. 


m 


2t 


D{s)ds 


m t 


(44) 


D{t) is then the time average of D{s) in the interval (0, t). 
The inconvenience of the use of D{t) is discussed in next 
section. Fig. ©• 

For completeness, note that, with a simple change of 
variables the z-GFPE, Eq. (IMll . can be written in the 
usual form of a diffusion equation. Defining the relative 
position Z{t) = z{t) —{z{t))^^ = z{t) — zo — (l)z{t) and using 
the relationship 


/dp{Z, t\Zo)\ 

fdp{z,t\zo)\ 

V dt 2 z{t) 

\ dt ) 


dt) 


fdp{z,t\zo)\ 


ip{t) = ipy{t) + 4>v{t), 


(49) 


in which, the random velocity term is given by a colored 
noise function Lpv{t) satisfying {(py{t))ii = 0 and the ex¬ 
ternal force was assumed to be time dependent but, con¬ 
stant in space. Therefore the moments of the ivspace 
probability density p{v,t \ vq) should be closely related 
to those of p{z, t I Zo,vo). In fact, using Chandrasekhar’s 
argument [l], [iSj , we show in Appendix A that p(u,t | vq) 
can be written as a normal Gaussian distribution, see Eq. 
(|A-6I) . Gonsequently, it satisfies a generalized master equa¬ 
tion identical to the velocity Rayleigh equation for the LE, 


2 dp{v,t I vo) \ 

V dt i 


dp{v,t I uo) 

- 

n I Vq) 


(50) 


but with the first and second moments generalized as 


= voXv(t) + (fvit), 


(51) 


= {plit)) 


KT 



4>l{t). 


(52) 


We shall refer to Eq. (1501) . together with (IMl) and (15^ as 
the generalized velocity Rayleigh-Fokker-Planck equation 
(v-GRFPE). 

With the definition of the total response function ip, Eq. 
dm), it is straightforward to obtain the relationships 

*(i) = -IJ{t)v{t) + P{t)>fi{t)+p{t), 

v(t) = -/3(t)v(t) + P{t)(j)y{t) + (j)^(t), 


(45) 
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(53) 

(54) 



















where 


m 


- ’Pyjt) 

v{t) - ’ 

-Xv{t)/Xv{t)- 


(55) 


Substituting this in the v-GRFPE, Eq. we then 

get the standard form of the velocity generalized Fokker- 
Planck equation (v- GFPE)[il,[il| 


dp{v,t I Up) 
dt 


+ 


I Vo)] 

nA.,^ d‘^p{v,t I vq) 

" ^ ’ dv^ 


(56) 


In which the drift and diffusion terms are 


^'drift(t) = V{t) + {Av)p, 


(57) 


and 

D^{t) = D^{t)+P{t)al{t) = + (AZ7),, (58) 

m 


where the shift in the drift and dispersion coefficients due 
to the presence of the time dependent external force can 
be written in Adelman’s notation [l^ as: 


(Au)^ 


u\ ^ r ^'>'(^) 


dt 


(59) 

(60) 


the velocity distribution is stationary, since the average of 
v^{t) over the probability Pi['c(t)] is 

{vHt)). = ^, (64) 

m 

This is consistent with the physical equipartition initial 
condition (uq)„ = k-gT/m, described in Appendix A, 

Ea. dAGB . 

3.3. Alternative view of the SDE 
Equation (l22ll is formally associated to the SDE 

dz{t) = v{t) dt + ify (t) dt. (65) 

However, since Eq. (l2^ is a well defined master equation, 
by inspection we can reinterpret it as resulting from the 
direct application of Ito’s lemmaj^ to the following SDE: 

dz{t) = v{t) dt + ^j2Dq{t) f{t) dt. (66) 

It has the same form of the drift term, but now the ran¬ 
dom contribution is given in terms of the standard Gaus¬ 
sian (5-correlated white noise f{t). It is weighed by a 
time dependent diffusion coefficient Dq(t), defined by Eq. 
(1^ . This result justihes the use of Hanggi’s type SDE, 
Eq. ra discussed in Appendix A, with f{t) = z{t), 
a{t) = v{t) = voXvit) + (fvit), bit) = Dq{t) and, Eq. (|A-5() 
with p{C,t) = p{z,t I zo,vo). This argument is also true 
for Eq. (1551) . which results from 


In the case of a free diffusing particle, (j)v{t) = 0, and 
this reduces to the standard form of the Fokker-Planck 
equation, Eq. (IA-81) . as derived by Adelman|13 A similar 
result was suggested by Hanggi and Talkner[9| but they 
omitted the last diffusive term, since their FDT did not 
contain the extra force term pointed out in Appendix A, 
Eq. (IA7M . 

The velocity probability density pi[v{t)] = p(v,t), irre¬ 
spective of the initial velocity, is obtained averaging over 
Pi[uo], the distribution of the vq, i.e 

pOC 

Pi[v{t)]= p{v,t \ vo)pi[vo]dvo. (61) 


Using the Gaussian distributions with the proper standard 
deviations, Eq. (IA-61) ioi p{v,t \ uq), and Ea. dlini) forpi[no], 
we obtain 


Pi[v{t)] = 




: exp 


[u(t) - 


2(7?, 


(62) 


2 2 I 2u\ 

= <^u+ -t), 

m 

= ^-cflit). (63) 

m 

Since fiyit) = 0 at t = 0, crJ(O) reduces to k^T/m , as 
expected. However, for t > 0 the pi[z;(t)] has a drifting 
term due to the effect of the external force. Nevertheless, 


dz{t) = (j)y{t)dt + \j2Dz{t) ^{t) dt. (67) 

This result corresponds to the use of Eq. with 

f(t) = z{t), a{t) = <l)v{t), h{t) = D^{t), and Eq. (lA-SP 
withp(C,t) =p{z,t I 2:0). 

Even though the dynamics of the SDEs, Eqs. (1551) . (|55)) 
and dnzi), are different, they are equivalent in the sense 
that their statistical properties should be identical. In fact, 
from Eq. dMl) 

I7(s) ds + 

Jo 

and recalling that {^{t)) = 0 and {£,{t)^{s)) = 6{t — s), we 
have, after taking the f{s) distribution average 

(^(t))^ = Zq + v{s) ds = (^z{t)'j{69) 

(z(t)2 - z(t)2)^ = 2^ Dq{s)ds = al{t), (70) 

in agreement with Eqs. (EH) and (EH)- Therefore, equa¬ 
tion (1551) also describes the position of a 1-D Brownian 
particle evolving with time in the GLE picture. This new 
SDE is simpler than the original, having the form of the 
Smoluchowski equation, but is valid for all values of the 
static friction coefficient and not limited to the condition 
imposed by the high friction limit approximation. It re¬ 
duces to the static LE for large A. In the same manner. 


^2Dqis)as)ds, ( 68 ) 


z{t) -Zq = [ 

Jo 
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Eq. (113 is in agreement with Eqs. OT and dll]) with 
{z{t))^^ — zo = Hence, Eqs. (155)) and (1571) are 

alternative Langevin processes that describe the position 
z(t) of a non-markovian particle with probability densities 
p{z,t I zo,vo) and p{z,t \ zq) respectively, related to the 
colored noise process, Eo. dhhl) . Their simplicity could be 
useful in the description of more complex systems. We 
point out that the velocity distributions associated to the 
above three SDE must be different in order to be consistent 
with the scaled ^ space. 


4. Analytical solution of the position z-GRFPE 
and z-GFPE for a fluid under a constant force. 


To describe a real viscoelastic fluid one needs to specify 
its bulk long time self diffusion constant Dq, besides other 
molecular parameters as the mass m and the temperature 
T. The bulk static friction coefficient is then given by the 
Sutherland-Einstein relationship 7 = k^T/D q, which is 
equivalent to the definition D^a = k^T/m with a = 7 /m. 
To solve the GLE, the friction kernel T{t) needs to be 
specified. As a first order approximation, we consider the 
standard exponential decaying kernel 

r(t)=aAe“^*. (71) 


This is convenient because it depends only on two 
parameters, and renders simple analytical relationships. 
Here, the frequency a measures the magnitude of the static 
friction, while 1/A is the memory relaxation timej^, [ 2 ^ . 
For this approximation, the characteristic Green function 
Xv{t) and the auxiliary functions Xzit)^ and x(t), Eqs. ([5]), 
(dlHl) respectively, can be evaluated analytically by solving 
the inverse Laplace transform of r(fc)[Ilj, 


Xvit) 

aXz{t) 

ax{t) 


= e -^‘/2 


cos (uj t) + -— sin (oj t) 
^ ' 2uj ^ ' 


= t 4 


1 - e -^‘/2 

X — a 


a A 
A —3 a 


/ \ A—2a . , . 

cos (w t) H— -sm (tu t) 

2 LJ 


e ^‘/^cos (wt) — 1 


2 aw 


e sin (wt). 


(72) 

,(73) 

(74) 


where w is an auxiliary frequency, 2 w = -^/A (4 a — A). 

For the case of a constant external force /)z(t) = Op x(f) 
and 4>v{t) = apX^(t), where Op = v^a = F/m measures 
the magnitude of the force. The parameters for the prob¬ 
ability density p{z,t \ Zo,vo) in the z-GRFPE are 

z{t) = zo + v„Xzit) + a^x{t), (75) 

a^{t) = Doa[2x{t) - xlit)] - alx^it), (76) 

Dq{t) = DQaXz{t)[l-Xv{t)]-alxz{t)x{i)- (77) 


While the parameters for p{z,t \ zq) in the z-GFPE are 
obtained as 


(^( 0 ).o = ^o+apX(0, 


(78) 


al{t) = 2Doax{t) - (79) 

-Dz(i) = Daaxzit) - alxz{t)x{t). (80) 

A central result is that the MSD parameters are indepen¬ 
dent of the applied external force 


a^{t) = 

2Doax{t), 

(81) 

m = 

Do ^ , 

(82) 

m = 

Dq — [cos u}t+— —^ sinwf]. 

(83) 

2lu 

The moments 

for the velocity v-GRFPE probability 

den- 

sity p{v,t 1 Vo) 

are 


v{t) = 

voXvit) +apXzit), 

(84) 

= 

[1 xl{t)] alxl{t). 

m j 

(85) 

Du{t) = 

-DoaXvit)Xvit) - alit) Xv{t)Xz{t), 

( 86 ) 

while those for 

Pi[v, t\ they should be 


(^(i))„o = 

= a^Xz{t), 

(87) 



( 88 ) 



f/ps 


Figure 1: The free diffusion reduced position and velocity fluctuation 
coefficients D*{t) = D{t)/DQ and TJj(i) = Du{t)/{DQa^) versus 
time t for an argon-like fluid. Solid lines are the GLE prediction for 
several A values. Dashed curves correspond to LE and superimpose 
to the GLE curve for A > 50ps“^. Also shown is the apparent D{t) 
function for A = 3ps“^ (dotted line). 

In Fig. © we compare the behavior of the fluctua¬ 
tion coefficients for position, D{t), Eq. (155)) . and velocity, 
Du{t), Eq. (1851) . as a function of time. The figure is for 
Op = 0 , but finite values of the force only affects the val¬ 
ues of the velocity U„(t). The GLE results are shown for 
several values of A. Particularly, for a very large A > 50, 
the curves match the classical LE (dashed curves), given 
by Eqs. lESl) and (IB-811 of Appendix B. For smaller A, 
oscillations are present and important discrepancies with 
LE are observed at times of the order of 1/a. As f ap¬ 
proaches zero, the ballistic behavior is obtained. We can 
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notice that the time scale of the velocity is much shorter 
than that of the position. It is an indication that the ve¬ 
locity field attains an earlier canonical distribution. Also 
shown is the apparent D{t) function for the low A = 3 ps“^ 
(dotted line). We can see that D{t) depicts only the time 
average of the diffusivity D{t) for the same A. At low times 
it misses the important oscillations present in the actual 
D{t). It requires extremely large times to converge to the 
correct limit Dq. Therefore its use is misleading. 


5. GLE Mean First Passage and Persistence Time. 


So far we have dealt with unbounded particles. In this 
section we shall take a look at the diffusion of constrained 
particles. First passage time distributions have been ex¬ 
tensively studied due to its related applications to diffusion 
controlled kinetics and biological processes 3ll-34|. Start¬ 
ing from the high friction limit diffusion equation, many 
physically valid boundary conditions have been applied, 
leading to well known analytical results. Based on the 
simple representation obtained in previous section, here 
we shall develop the GLE theory of first passage times of 
constrained viscoelastic fluid particles, under an external 
force. So, the z-GFPE, Eq. (I35L must be solved with the 
appropriate boundary conditions for one and two absorb¬ 
ing bounding barriers. 


5.1. Diffusion Next to a Single Absorbing Barrier 

Let us first consider a fluid particle initially at zq, diffus¬ 
ing under a force F'(t), but constrained in the z direction 
by an absorbing barrier located at position b > zq. The 
absorbing property of the barrier requires that the bound¬ 
ary condition p{b, t \ zq) = 0 is satisfied. This can be 
fulfilled representing the bounded probability by the stan¬ 
dard combination of unbounded p^{zff \ Zq) functions [i|, 
namely 

p{z, t\ Zq) = p°{z, t\ Zq)- p° {2b - Z,t\ Zg) . (89) 


Integration over z gives, after some algebra, the survival 
probability 


G{zo,t) 


p{z,t I Zq) dz, 


erf 


Zbjt) 1 


(90) 


where Zb{t) = b— Zq — 4>z{t), and erf[z] is the error function 
of z. 

Defining /i(t, Zo)=—dG{t, zff^jdt as the first passage time 
distribution, we get 


h{t,zo) 


X 


2[Zb{t)Dfft) + (j)y{t)al{t)] 


exp 


(t) 

zm 1 

2al{t). ■ 


(91) 


In the absence of an external force, this reduces to 

<^Xz{t){b - Zq) 


[h{t,Zo)]F=0 = 


ffA-irDo [axp(t)]^'^^ 
{b- zo)'^ ^ 


exp 


4DoaXp(t) 


(92) 


Equation (lOTl) is an exact closed relationship for the dy¬ 
namic generalized Langevin equation of a viscous fluid un¬ 
der a time-dependent space-fixed finite external force. For 
a time-constant force and an exponential kernel, it is an¬ 
alytical, since the functions Xv {t), Xf (^) i 4'v (t), cr^ {t ), and 
Dz{t) were given analytically in previous sections. Partic¬ 
ularly, in the high friction or Smoluchowski limit (A —> oo, 
a —>■ oo, and Dfft) —>■ Dq), Eq. (IMl) becomes for a con¬ 
stant force 


[^(^:^o)]se = 


X 


h- zo + vff 
i/47riAo 


exp 


{b - zq - 
4iAot 


(93) 


where, as we defined before, v^t = FDff/k^T. Note that 
Eq. (IMll differs slightly from the Smoluchowski result for 
a constant force obtained by Hu et aL[^ where the drift 
contribution is missing in the numerator. Additionally, 
our contribution contains a new term, quadratic in the 
force. 


5.2. Diffusion Between Two Absorbing Barriers 

Another physical situation widely studied in the SE limit 
corresponds to the self -diffusion of a viscous fluid within 
an absorbing virtual slab of molecular dimensions. Let z = 
0 and z = L be the position of two absorbing boundaries 
(A-A) of a virtual slab of length L. The absorbing bound¬ 
ary conditions require that p{0,t \ zq) = p{L,t \ zq) = 0. 
For a time dependent drift velocity fiv{t), due to the exter¬ 
nal force F{t), the solution of Eq. (1551) . can be obtained by 
separation of variables UM- Using the initial condition 
p{z, t = 0 I Zq) = 5{z — Zq) , we obtain 


O oo 

p{z,t\zo) = - 2^sm(—z)sm(—Z q) 


X exp 


(94) 


where we introduced the dimensionless time dependent 
Peclet number. 


Pe(t) 


fiv{t)L 

2Dffty 

Pe° 

l-{?ZffD,ax{t) 


(95) 

(96) 


which measures the importance of the mass transfer due to 
the external drifting force, relative to that due to diffusion. 



























Here Pe° = Pe(< = 0) = vfL/2Dq. The function r„(t) is 
a re-scaled time given by 


Tn{t) 


[n^ + Pe^ (t)] f* 

Pe"(t) Jo 

7 r^-|-Pe^(t)] 

4i?oPe"(t) 


ds, 

2(Pe°)V2(t) 


In 


1 - 


L 2 


(97) 

(98) 


The survival probability for a particle to remain in this 
region having started to diffuse at zq > 0 is given by 


Gizo,t) = / p{z,t I zo)dz. 
Jo 

A direct integration gives 


(99) 


n=l 


X exp 


n2 TfZ _|_ Pe2(t) 
=0 


- —^0 - 


( 100 ) 


Note that using the definition of G{zo,t) in the backward 
version of Eq. dMl), we find that the survival probability 
G{zo,t) must obey the differential equation 


dG{zo,t) ^ dG{zo,t) d^G[zo,t) 

G(0,t) = G(T,t) = 0, 

G(zo,0) = 1, (101) 

whose solution is just Eq. (HOOD . 

The distribution of the first passage times, h{zo,t) can 
be evaluated using the fact that it is equal to the net flux 
reaching the absorbing boundaries 


h{zo,t) = - 


dG{zo,t) 


dt 


= D, 


LV oz Jl \ oz Joi 


then 


h{zo,t) = 2^ 


n=l 

X sin 


n:D^{t) 

““L 2 


/nTT \ 


exp 


l_eP«d)(-l) 
-^^0-^Pn(4 (103) 


The persistence probability P{t) that the particle is still 
in the slab ( 0 ,L) at time t, irrespective of its initial po¬ 
sition, is obtained averaging the survival probability over 
the distribution of initial positions 


P{t) = {G{zo,t))^„= G{zQ,t) g{zQ) dzo, (104) 
0 

where g{zo) is the one-particle radial distribution function 
of the fluid, normalized in the slab [0,L]. In the virtual 
layer model simulations [ 3 , ^^e external mean force is 


assumed to be a constant within the slab, so the poten¬ 
tial of mean force is W{zo) = —Fzq and in terms of the 
parameter Pe° 


gi^o) = -pNgexp 


2 Pe'' 

CT 


-zo 


(105) 


With Ng = 2Pe°/[exp (2Pe°) — l] as the normalization 
factor. Integrating over zq, the persistence probability in 
the absorbing slab is 


P{t) = 2Ngj2 


2 2 
n TT 


I-(-I)"e 


Pe“ 


n—1 


r2 7r2 -p Pe^^ 


I-(-I)"eP®(*) 

n2 7r2 _j_ Pe^(t) 


exp 


^2 Pn(^) 


(106) 


The mean first passa ge t ime tMPp('2o) (MFPT) is the first 
moment of h{zo, t)[i,Mn 


nOO nOO 

^mpp(2^o)= / th{zo,t)dt= / G{zo,t)dt. (107) 
Jo Jo 


Its average over the distribution of the initial positions is 
called the persistence ti? 7 ie[ 3 ^, denoted by r. For absorb¬ 
ing boundaries it is also referred to as the mean exit time. 
In terms of P(t), t is given by 


T = {tMPp{zo))zo = p{t)dt. (108) 

Jo 

That is, the persistence time is just the normalization con¬ 
stant of the persistence probability. 

The mean square displacement for particles persisting in 
the absorbing slab at time t, (r\{t), is a bounded function. 
Formally, it should be obtained integrating {z{t) — zoY over 
z(t) and Zo with a bounded probability distribution given 
by Eq. ([Ml). However, another simpler expression can be 
computed in terms of the unbounded diffusion coefficient 
and the survival probability P(t) 



D{s)P{s)ds. 


(109) 


5.3. A-A Diffusion for Small Pe Numbers 

Peclet numbers, so defined above, are in general nonlin¬ 
ear functions of time, which depend on the functional form 
of the external force F{t). Therefore, the time integrations 
needed to evaluate Impp ( 2 ^ 0)1 its average r and <j\{t), have 
to be carried out numerically according to the exact GLE 
closed relationships given by Eqs. (I107D . (I108p . and (I109|) . 
for a given value of A. 

To obtain analytical results, let us consider the limit 
of small Peclet numbers. For a constant and relatively 
small force, the Peclet function is a constant to quadratic 
order in the force and we can approximate Pe(t) ~ Pe°. 
Denoting it simply as Pe, we can simplify 

Tn{t) = TT^ -f Pe^)ax(t), 
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( 110 ) 











































so we can readily evaluate 


al{t) = m^L 


n—1 


7r2[l-(-l)’"eP®]^ 

[n^ TT^ + Pe^]^ 


1 — exp 


Pe" 


2L2 


a^{t) 


For a bulk fluid with zero force, Pe = 0, 
Tnit) = n^TT^axit), 


( 111 ) 


( 112 ) 


and 


n—1 


^ exp [ - ^^Ij^Doaxit)]. 


(113) 


This is consistent with the persistence probability obtained 
from Eq, (110611 for Pe = 0 




n—1 


TT^ 

L2 


Doaxit) .(114) 


Since the particles in the slab are continuously being ab¬ 
sorbed at the boundary, P{t) decays fast to zero and the 
(j\{t) reaches a plateau 1 ?The last is a geometrical 
factor associated to the phenomenological description in¬ 
herent to the GLE approach. 

To the best of our knowledge, all attempts to study a 
fluid diffusing within absorbing barriers are based on the 
Smoluchowski or high friction approximation [i , ^ ^ 37, 
[sil which can be obtained as a particular case of our gen¬ 
eral equations. As A —oo we get the static Langevin 
approximation, using q:Xle( 0 = (at — 1 -I- e““*)/a, 
while as at ^ 1 we get the Smoluchowski limit, using 
aXssit) = t and Dz{t) = D{t) = Dq. In fact, for small 
constant Pe numbers, dropping quadratic terms, we let 
Pe = FL/2mDQ, and we can readily integrate G{zo,t) 
over time to get the following analytical expressions in the 
SE limiting approximation 


.SE 7 . ^ 

MPP^'^Oj jj / ^ 

^ n—1 


TT [1 - (-l)”eP®] 
[71,2 7r2 -I- Pe2]2 




X sin(n7r—) exp[—Pe—], 

Lj 1j 

AL^Ng ^n2 7r2[l-(-l)"eP®] 


Dn 


a=i [ 7r2-I-Pe^] 


21 3 


(115) 


(116) 


[ai(t)]SE = 2Dot^^ - AL^NgY, 


n2 7r2[l-(-l)"eP®]' 


exp - 


„^i [ n2 7r2 -H Pe^] ^ 

n? 7r2 -I- Pe^ 


L2 


-Dot] , 


(117) 


Eq. (11151) is very close to the expression obtained by 
Gitterman [32| for the MFPT, which nevertheless misses 
the finite exponential terms containing Pe. 


For zero force, the sums appearing in Eqs. dm]) 
and (I116|) become the Riemann type, namely, 
[1 - (-l)"]sin(n7r2:o/A)/n^ = - zo)zo/AL‘^ 

and [1 — “ '^^/48, so we get the well 

known results [321 [SSj, l39|| 


G^^{zo,t) 


1 — (—1)" . /nTT 


J F=0 


= 2^^-^^ sm 

' n TT 


n=l 


fniT \ 


X exp 


-\ — ]Dot 


(118) 




n—1 




[Cp(^0)]f=0 = 


X exp 
(L - zo)zo 


-[-)Dot 


2Do 


^SEi 


12iAn’ 


(119) 

( 120 ) 
( 121 ) 


Kilto = 

n—1 


X exp 


Do t 


( 122 ) 



^/ps 


Figure 2: The time dependent cr‘^(t) = MSDij{t) in an absorbing 
slab of length L = 1.66 A, located in the bulk (a^ = 0), normalized 
by its saturation value. GLE (solid line) for A = 3 ps“^, a = 5.5 ps“^ 
in Eq. Ill 1311 . LE (dashed line) and SE Wot-dashed line). The black 
dots correspond to the MD simulation^. 


To test the above equations, we used the MD simulation 
data for LJ-argon fully described elsewhere!^. In Fig. (I2|) 
we show the time dependent MSD for particles constrained 
in an absorbing slab of length L = 1.66 A, located in the 
bulk (op = 0). Values are normalized by its saturation 
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value, namely MSD\(t) = a^it)/(7\{oo). The theoreti¬ 
cal predictions are compared to the MD simulation (black 
dots)0. An almost perfect match was obtained in the 
GLE (solid line) with an exponential kernel, according to 
Eqs. (Illdl) and (I7H) . for {A = 3ps“^, a = 5.5ps“^}. The 
LE (dashed line) approximation predicts the initial bal¬ 
listic behavior but fails to give the characteristic s shape 
of the saturation MD curve. The SE (dot-dashed line) 
approximation predicts only roughly the qualitative be¬ 
havior. 

6. Conclusions. 

We have derived analytical formulas for the solution of 
the non static or generalized Langevin equation (GLE) 
of a particle in a viscoelastic fluid, under the action of 
a time dependent external field. Emphasis was done in 
unifying the notation and existing different forms of the 
master equations, not just for the velocity distribution 
P{v,t I woiO) and its dispersion coefficient Du{t), well 
discussed in text books 0, [13 ; but also for the condi¬ 
tional position probability densities P{z,t \ zoi'i'OiO) and 
P{z,t I ^OjO) and their corresponding diffusivity coeffi¬ 
cients Dq{t), Dz{t) and the time dependent diffusion coef¬ 
ficient, D{t), associated to the MSD. In doing so, we have 
extended the results of Dufty[Il| and Rodriguez et al. [l3 
for an unbounded particle to include a time dependent 
external force. The GLE dynamics of bounded fluid par¬ 
ticles constrained by absorbing barriers was then obtained 
in terms of the results for the unbounded particles. 

We wrote all dynamical properties in terms of the funda¬ 
mental Green functions of the GLE x(t), Xz{t) and x„(t), 
which are the natural susceptibilities of the system and 
extended the fluctuation-dissipation theorem to include 
the external force. Since there is no solution for a gen¬ 
eral external field, we restricted the equations for the sim¬ 
ple constant force case. Thus, we found that for an ex¬ 
ponential decaying friction kernel the analytical formu¬ 
las obtained, like Eq. (nm, fit well the available exact 
molecular dynamics simulations. We tested our equations 
against molecular dynamics data of Argon atoms diffus¬ 
ing between two impenetrable walls separated 40 atomic 
diameters each other Q. In particular, at the center of 
the system, where the fluid behaves field free, we found 
that the mean square displacement in a slab of length 1.66 
A, normalized to its saturation value, almost matches the 
simulation. 

We also would like to point out that subsequently from 
this formal GFPE, we have obtained statistically equiva¬ 
lent stochastic differential equations for the position of the 
non-static motion, Eqs. (1551) and dnzi), where the noise 
is white, as in the static Langevin case, but with time 
dependent diffusion terms Dq{t) and Dz{t). The GFPE 
associated to these alternative non-static SDE have the 
properties expected for a Brownian particle. Namely, its 
stationary solution for the position probability densities 
p{z, t\zQ, vq, 0) and p{z, t\zo, 0) are Gaussian. In summary, 


it rigorously treat the Brownian motion in the configu¬ 
ration space, discarding the instantaneous velocity of the 
particle without resorting to a limit operation of any kind. 
The new SDE does not impose any restriction on the mag¬ 
nitude of the fluid friction coefficient. For the field free 
case, it has the same average and standard deviation as 
those prescribed by Ghandrasekhar[l3. Moreover, the 
Smoluchowski limit is guaranteed in this new approach. 
This limit was specially useful in problems where the fric¬ 
tion term is dominant, such as in the analysis of the dif¬ 
fusion coefficient in non-homogeneous fluids. They in¬ 
clude argon atoms restricted to diffuse between the force 
of two impenetrable walls [s^, biological systems [sij and 
ionic fluids 3^, 4^ . Our results could be applied to such 
systems with the interpretation that {t) is now a gener¬ 
alized non-static diffusion coefficient. That is, the existing 
approaches based on the Smoluchowski equation would be 
expanded to the domain of low friction coefficient. 

Our GLE describes diffusion that does not obey the nor¬ 
mal Brownian relationships. For instance, the apparent 
anomalous behavior shown by Eq. (82) is mesoscopically 
due to a retardation in the friction under a colored noise. 

Besides the system studied here, this method of analy¬ 
sis can be perfectly applied to exit time problems (3, con¬ 
trolled diffusion reactions [siils^ - diffusion between reflec¬ 
tive barriers [s^l, stochastic thermodynamics [dj] and feed 
back control in small systems (d^. to mention just a few. 
All of them have been analyzed from the perspective of 
Smoluchowski overdamped motion, so our procedure can 
be used to extend them to the demanding underdamped 
regime. 
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Appendix A: The GLE velocity fluctuation coeffi¬ 
cients. 

In this Appendix we give a unified view of the existent 
treatments of the velocity as a stochastic variable. This 
will led to the evaluation of the fluctuation coefficients 
associated with the GLE. 

Most of the mathematical treatments for a general 
stochastic variable or process ({t), have focused on a gen¬ 
eralized SDEs of the kind 

C{t) = a{C,t) + (A-I) 

where the drift a((, t) and the random force intensity 6(C, t) 
are explicit functions of time and the variable The exis¬ 
tence of formal master equations for this problem has been 
proven Q. The result, usually obtained via a functional 
Taylor series of its cumulants, shows that the process C,{t) 


II 





is markovian and the resulting master equation is Fokker- 
Planck-like in the Stratonovich sense i.e., 


dpiC,t) 

dt 




1 

2 5C 


KC,t) ■^KC,t)piC,t) 


(A-2) 


But we should point out that, using the relationship be¬ 
tween the Stratonovich drift velocity first moment and 
Ito first moment (^)j 


{Os = {a + l§^b\c,t), 

the corresponding FPE in the Ito sense is 


(A-3) 


dpiC,t) 

dt 




+ 


1 

2 ^ 


[b^{C,t)p{C,t)] ■ 


(A-4) 


When ({t) is the velocity v{t), Eq. (IA-411 is just a Rayleigh- 
type equation for the velocity probability p{v, t \ vq), given 
an initial velocity vq. In fact, for a{v,t) = a(t) = v{t) and 
h{v,t) = h{t) = r](t)/^(t) = yj2Du{t), independent of v{t), 
we get the master equation associated to the LE, Eq. o 
in the absence of an external field, i.e. F(z, t) = 0, namely 


dp{v,t I Vq) 

dt 


-f,. dp{v,t I z;o) 

- 

d‘^p{v,t I Vo) 


Du{t) 


dv"^ 


(A-5) 


where v{t) = (n(t)){ is the average of v{t) over the noise 
We have used the subscript u, instead of in the fluc¬ 
tuation coefficient D„, anticipating that the resulting Eq. 

has the standard form of the Kramers-Moyal master 
equation expansion, with a first drift parameter in terms 
of the rate of the first moment v and a second disper¬ 
sion parameter in terms of the rate of the second moment 
of u{t) = [v{t) — u], namely, 2Du{t) = du^{t)/dt. The 
appropriate application of Eq. (E3 when the stochastic 
variable C,{t) is set as the particle position z(t), is discussed 
in the text. 

Using the standard Chandrasekhar’s arguments [H [I^, 
the quantity u{t) must have the same Gaussian statistics 
as the noise, so the velocity conditional probability density 
which satisfies Eq. can be readily written as 


p{v,t I Vo) = 


\/2 7r(T„2(t) 


exp 


2o-2(t) 


(A-6) 


where da^{t)/dt = 2Du{t). Writing the Gaussian in the 
alternative form 


dp{v,t I Vq) 
dv{t) 


[u(t) - n(t)] 


p{v,t\vo), (A-7) 


we can rewrite Eq. (IMJ in Adelman’s Fokker-Planck form 

n 


dp{v,t I Vo) 
dt 


= m 


d[v(t)p{v,t I Vo)] 


Dtit) 


dv 

d‘^p{v,t I Vq) 

dv"^ 


(A-8) 


where /?(t) = —v{t)/v{t) and the apparent coefficient 
D^{t) is 

D^it) = Du{t) + 2 m [ Duis)ds. (A-9) 

JQ 

This form of the master equation for the generalized GLE, 
Eq.Q, is commonly referred to as the velocity space 
Fokker-Planck equation (v- FPE)[i,[Il,[l|. Next, we eval¬ 
uate the properties of the colored noise functions and the 
jump moments corresponding to the probability densities 
p(v,t\vo) and p{v,t) for the velocities in Chandrasekhar’s 
context [18|. 

According to Eqs. (HUD and (USD, for a Gaussian color 
noise the definition of the auxiliary velocity function 


u{t) = v{t) - v{t) = (fivit), 


(A-10) 


with v{t) = voXvit) + 4>v{t), ensures, a diffusion-like 
equation for the velocity conditional probability density 
p(u,t\vo) 


( 


dp{u,t\vo)\ _ ^ p(u,t\vo) 

dt )u{t) dv? 


(A-11) 


Equations (IA-411 . (IA-511 . (IA-811 . and (lA-llIl are all equiv¬ 
alent master equations in v-space for the GLE p(v,t \ vo)^ 
when the proper definitions of the diffusion coefficients 
Du{t) and D^{t) are used. The solution of Eq. (lA-lip 
in an unbounded space is the standard Gaussian 


p{u,t\vo) = 


1 




u^t) 


(A-12) 


The first moment is simply u(t) = {u(t))^ = 0 and 
(■^(O)r = = VoXvit)d-fj^vit)- The velocity drift <j)v{t) is 

defined in Eq. (USD and reduces to a^Xzit) for a constant 
force. 

The second moment, cr^it) = {u‘^{t))j^, can be obtained 
from Eq. (IAAOI) 


^lit) = {u\t)), = {plit)),. 


(A-13) 


The first identity, together with the solution of the GLE, 
Eqs. (USD and da, gives 

(^lit) = {v^it))^ - Voxlit) - (l^lit) - 2voMt)- (A-14) 

Since cr'^it) is independent of the initial conditions, one 
can average this identity over vo and use the physical sta¬ 
tionary and equipartition conditions 


((vHt)}^L = (^olo = 


KT 


m 


(A-15) 
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with (vq)^^ = 0, in order to get 

= — [1-xSW] - (A-16) 

m 

and therefore the diffusion like velocity coefficient Du{t) 
in Eq. (|A-llll is = (l/2)dcr2(<)/(it 

Du{t) = (A-17) 

m 

where 

(A-18) 



Figure 3: The free diffusion two-time correlation function C(t, s) = 
{t)p^ (s)) vs (t — s), scaled by {Dq a). Graphs (a) and (b) are for 
a null external force at two initial times: (a) a low s = 0.09 ps and 
(b) a higher s = 0.25ps. Graphs (b), (c) and (d) show the effect of 
a constant external force: a-^ in A/ps^ of 0, 1.5 and 3.0 respectively. 
In all graphs a low and a high values of A were used for GLE (solid 
lines), 3ps“^ (labeled) and 25ps~^. Dashed curves correspond to 
LE. Q = 5.5 ps“^. 


These results can be obtained also from the second iden¬ 
tity of Eq. (|A-13p involving the color noise force function 
Pv{t) defined in Eq. (fT^ 11, 1^ 


Xv{t-t') dt'f Xv{t-t"){R{t')R{t''))dt''.{A-19) 
Jo Jo 

The Gaussian R-noise two time correlation appearing 
in this expression must satisfy a physical consistency rela¬ 
tion, which is widely known as the fluctuation-dissipation 
theorem (FDT). For a time dependent external force F{t) 

+ F{t')F{t")'^ = 

-t" \), (A-20) 

m 

k^TmT{\ t'-t" I) 

F{t')F{t"). (A-21) 


—,{{R(t’)Rr)), 


This relationship is derived from Eqs. @ and ©HO], 
through the correlation of the dissipative force r(| t' — t" \ 
) with both the internal and external fluctuation force 
Rit) + F{t), assuming the stationary and equipartition 
conditions, therefore it is equivalent to Ea. (IA-15|) . Here 
we have an extra drift term due to the external force which 
is often neglected in the literature 4^, M\- This theorem 
is useful for evaluating the two-time correlation function 
C{t, s) = ; it is given by: 

C{t,s) = / Xv{t-t')dt' 

^ Jo 

X ^ xv[s - t") (^R{t')R{t")^ dt". (A-22) 
It was nicely shown by Adelman[l^ and Fox[I3 that 

fx,{t - 1 ') dt' / x„(s - 1 ") r(| t' - 1 " I) dt" 

Jo Jo 

= Xvi\t-s \) - XviOXvis)- (A-23) 


Using this, together with Eq. (IA-20|1 . in Ea. (IA-19l) one 
obtains the general result 


C{t,s) 


— [x«(U - s I) - X« WXi-(s)] 

m 

(l)y{t)4>v{s). (A-24) 


Using t = s and Xt,(0) = 1, Eq. (jA-241) reduces to Eq. 
(IA-161) as expected. 

In Fig. ©, the two time correlation function C{t,s) is 
shown as a function of the time difference (t — s), for two 
initial times s = 0.25 ps (top) and 1.5 ps (bottom) and two 
values of the frequency A as labeled. Static Langevin the¬ 
ory (LE) gives a single exponential decay (dotted curve), 
as derived in Appendix B, Eq. (IB-121) . For A as large 
as 25ps“^, GLE already approaches the LE. Since the 
frequency w was defined as a real positive quantity for 
4a > A, the susceptibility function XviO is an oscillatory 
function of time for low values of A. Hence a negative 
correlation is expected from the definition of C(t, s), Eq. 
(|A-24I) . particularly for the higher values of the external 
force. 


Appendix B: The static-LE and high friction-SE 
limits. 

In this section, we establish the connection of our gen¬ 
eral results with the commonly used approximations in the 
static and high friction limits. To find the drift and diffu¬ 
sion terms in the static LE for a constant external force, 
we apply the procedure of Sections and dH) to Eq. dl]). 
The result is similar to Eq.dH]), with A —?► oo. The char¬ 
acteristic functions become Xv^ = ay^^ = 1 — e““*, 
and a^x^^ = -|- at — 1. Furthermore, 

^le(^) =''^oe““‘ + ^p(l-e"“*), (B-1) 
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= Do a (1 - e-2“‘) - (b_2) 

[Duit)]^E = Doa^e“^“* - z;pa(l - e““*)e““*. (B-3) 

The moments associated to the probability density 
p{z,t I 2;o,r’o) are 

mr = z„ + ^(l-e-“‘) 

a 

+ ^(e-“* + at-l), (B-4) 

a 

= ^[2at-3 + 4e-“‘-e-2“‘] 

- ^(e-“‘+at-l)^ (B-5) 

For = 0 this reduces to the textbook result [I|. The z;o- 
average probability p^^(z,t\z^) given by Eq. (ISTl) . evolves 
in time as an equation similar to Eg. (1351) with fluctuation 
moments: 


mCg = 

■^0 + “ -h at - l) 

a 

(B-6) 

Wlit)]LE 

= -f af-l) 

a 

- 4(e-“‘ + af-l)/ 

(B-7) 

[D.{t)]EE 

= D„(l-e-“5 



X [1 («“*+“' 1)]. 

(B-8) 


The plain MSD{t) and diffusivity coefficient D{t) are 


9 Jl 

[MSD(f)U = ^(e-“‘+af-l), 

(B-9) 

[D(f)U =Do(l-e-“5. 

(B-10) 

The apparent diffusivity D(t) is 


[5(t)U = ^(e-“‘ + at-l), 

(B-11) 

corresponding to the colored noise (p^E (t), whose two-time 
correlation function for f > ,s is iust the limit of Eg. (IA-24II 

for A —?► oo 


- u2(l-e-“‘)(l-e-“"). (B-12) 


All of this agrees with the result of applying Eq. (TT^ 
consistently, with rLji;(t) = a d{t) and 

-^{R{t)R{s))EE = — (B-13) 

m 

respectively. We notice that both, the LE and its general¬ 
ization GLE, have a Gaussian probability distribution. In 


particular, for the static LE, Eq. (El]) with zero external 
force, i.e. = 0, matches Eq. (171) of Ghandrasekhar pa¬ 
per on Brownian motion of a 1-D free particle[l^. There¬ 
fore, at large times {t) — 2Dgt. The so called high fric¬ 
tion limit (HFL) corresponds to the limit a —>■ oo. Phys¬ 
ically, it is consistent with a small relaxation time 1/A of 
the friction retardation kernel. So, the HFL of the GLE 
requires that we take first the static LE limit, A —>■ oo, and 
a —> oo, afterwards. Since, the first limit corresponds to 
the static LE, then from Eqs. dEl) and (iBTa . the HFL 
gives Fhpl(^) = br ^hpl(^) = Rq- Therefore, Eq. (IHHl) 
becomes under this condition 

^=VE + V^om- (B-14) 

Likewise, Eq. ED) reduces to the so called Smolu- 
chowski equation (SE). The probability density PjjpL(z,t) 
is a Gaussian with the moments z(t)jjpL = Zg + Vpt and 
cr^pt (0 = 2 Dj, t. Another way to get Eq. (jB-141) is by 
taking first the limit A —>■ oo to Eq. (IA-221) and then 
the limit of large a. As expected, the hrst result repro¬ 
duces Eq. (IB-1211 , and then for large friction, the corre¬ 
lation {ip{t)ip{s)) is a Dirac delta function and the noise 
ip{t) transforms into white noise in this limitQ. Similarly, 
the drift reduces to Wp and Eq. (1551) becomes Eq. 

(iBHl) . 


References 

References 

[1] D. A. McQuarrie, Statistical Mechanics, Harper Row, New 
York, 1976. 

[2] J. E. Keizer, Statistical Thermodynamics of Nonequilibrium 
Processes, Springer-Verlag, New York, 1987. 

[3] L. E. Reichl, A Modern Course in Statistical Physics, 2nd Edi¬ 
tion, John Wiley & Sons, New York, 1998. 

[4] P. J. Colmenares, F. Lopez, W. Olivares-Rivas, Molecular Dy¬ 
namics and Analytical Langevin Equation Approach for the 
Self-Diffusion Constant in Anisotropic Fluids, Phys. Rev. E 80 
(2009) 061123. doi:10.1103/PhysRevE.80.061123 

[5] C. W. Gardiner, Handbook of Stochastic for Physics, Chemistry 
and the Natural Sciences, 2nd Edition, Springer-Verlag, Berlin 
Heidelberg New York Tokyo, 1985. 

[6] J. P. Boom, S. Yip, Molecular Hydrodynamics, Dover, New 
York, 1991. 

[7] H. Risken, The Fokker-Planck Equation. Method of Solutions 
and Applications, Springer-Verlag, Berlin Heidelberg New York 
Tokyo, 1984. 

[8] P. Hanggi, Correlation Functions and Masterequations of Gen¬ 
eralized (Non—Markovian) Langevin Equations, Z. Physik B 31 
(1978) 407-415. doi:10.1007/BF01351552 

[9] P. Hanggi, P. Talkner, On the Equivalence of Time- 
Gonvolutionless Master Equations and Generalized 
Langevin Equations, Phys. Lett. A 68 (1) (1978) 9-11. 
doi:10.1016/0375-9601(78)90740-5 

[10] T. S. Chow, J. J. Hermans, Effect of Inertia on the Brownian 
Motion of Rigid Particles in a Viscous Fluid, J. Chem. Phys. 56 
(1972) 3150-3154. doi:10.1063/1.1677653 

[11] J. W. Dufty, Gaussian Model for Fluctuation of a Brownian Par¬ 
ticle, Phys. Fluids 17 (1974) 328-333. doi: 10.1063/1.1694718 

[12] S. A. Adelman, Fokker-Planck Equations for Simple Non- 
Markovian Systems, J. Chem. Phys. 64 (1976) 124-130. 
doi:10.1063/1.431961 


14 


















[13] R. F. Fox, Gaussian Stochastic Processes in Physics, Phys. Rep. 
48 (3) (1978) 179-283. doi:10.1016/0370-1573(78)90145-X 

[14] V. S. Volkov, V. N. Pokrovsky, Generalized Fokker-Planck 
Equation for NonMarkovian Processes, J. Math. Phys. 24 (1983) 
267-270. doi:10.1063/1.525701 

[15] F. R. Rodriguez, E. Salinas-Rodriguez, Brownian Motion and 
Gorrelation Functions in a Viscoelastic Fluid, J. Phys. A: Math. 
Gen. 21 (1988) 2121-2130. doi:10.1088/0305-4470/21/9/025 

[16] A. A. Budini, M. O. Caceres, Functional Gharacterization of 
Generalized Langevin Equations, J. Phys. A: Math. Gen. 37 
(2004) 5959-5981. doi:10.1088/0305-4470/37/23/002 

[17] E. Darvea, J. Solomon, A. Kiab, Computing Generalized 
Langevin Equations and Generalized Fokker-Planck Equations, 
Proc. Natl. Acad. Sci. U.S.A. 106 (27) (2009) 10884-10889. 
doi:10.1073/pnas.0902633106 

[18] S. Chandrasekhar, Stochastic Problems in Physics 

and Astronomy, Rev. Mod. Phys. 15 (1) (1943) 1-89. 

doi:10.1103/RevModPhys.15.1 

[19] R. Kubo, Stochastic Liouville Equations, J. Math. Phys. 4 (2) 
(1963) 174-184. doi:10.1063/1.1703941 

[20] N. G. van Kampen, Stochastic Processes in Physics and Chem¬ 
istry, North-Holland, Amsterdan, 1981. 

[21] J. M. Sancho, M. S. Miguel, S. L. Katz, J. D. Gunton, Anal 3 rtical 
and Numerical Studies of Multiplicative Noise, Phys. Rev. A 
26 (3) (1982) 1589-1609. doi:10.1103/PhysRevA.26.1589 

[22] P. Liu, E. Harder, B. J. Berne, On the Calculation of Diffusion 
Coefficients in Confined Fluids and Interfaces with an Applica¬ 
tion to the Liquid-Vapor Interface of Water, J. Phys. Chem. B 
108 (2004) 6595-6602. doi:10.1063/1.2841128 

[23] P. J. Colmenares, Role of a Stochastic Friction Coefficient 
in Open Channel Noise, J. Theor. Biol. 161 (1993) 175-198. 
doi:10.1006/jtbi.1993.1049 

[24] V. S. Volkov, A. 1. Leonov, NonMarkovian Brownian Motion 
in a Viscoelastic Fluid, J. Chem. Phys. 104 (1996) 5922-5931. 
doi:10.1063/1.471324 

[25] N. G. van Kampen, Stochastic Differential Equations, Phys. 
Rep. 24 (1976) 171-228. doi:10.1016/0370-1573(76)90029-6 

[26] K. Furutzu, On the Statistical Theory of Electromagnetic 
Waves in a Fluctuating Medium (I), J. Res. Natl. Inst. Stand. 
Technol. 67D (3) (1963) 303-323. doi: 10.6028/jres.067D.034 

[27] E. A. Novikov, Functionals and the Random-Force Method in 
Turbulence Theory (Euler Velocity Field Described by Ran¬ 
dom Forces Method, Using Lagrangian Representation of Tur¬ 
bulence), Zh. Exper. Teor. Fiz. 47(5), 1919 (1964), Sov. Phys. 
JETP 20 (5) (1965) 1290. 

[28] M. D. Donsker, On Function Space Integrals, in: W. T. Mar¬ 
tin, 1. E. Segal (Eds.), Analysis in Function Space, MIT Press, 
Cambridge, Mass., 1964, pp. 17-30. 

[29] V. B. Gnedenko, The Theory of Probability, MIR Publishers, 
Moscow, 1978. 

[30] B. J. Berne, G. D. Harp, On the Calculation of Time Correla¬ 
tion Functions, in: I. Prigogine, S. A. Rice (Eds.), Advance in 
Chemical Physics, Vol. XVH, John Wiley &; Sons, New York, 
1970, pp. 64-224. doi:10.1002/9780470143636.ch3 

[31] A. Szabo, K. Schulten, Z. Schulten, First Passage Time Ap¬ 
proach to Diffusion Controlled Reactions, J. Chem. Phys. 72 (8) 
(1980) 4350-4357. doi:10.1063/1.439715 

[32] M. Gitterman, Mean First Passage Time for Anoma¬ 
lous Diffusion, Phys. Rev. E 62 (5) (2000) 6065-6070. 

doi:1063-651X/2000/62(5)/6065(6) 

[33] G. Hummer, A. Szabo, Kinetics from Nonequilibrium Single- 
Molecule Pulling Experiments, Biophys. J. 85 (2003) 5-15. 
doi:10.1016/S0006-3495(03)74449-X 

[34] Z. Hu, L. Cheng, B. J. Berne, First Passage Time Distribution in 
Stochastic Processes with Moving and Static Absorbing Bound¬ 
aries with Applicaton to Biological Rupture Experiments, J. 
Chem. Phys. 133 (2010) 034105. doi: 10.1063/1.3456556 

[35] P. J. Colmenares, W. Olivares-Rivas, Smoluchowski Hypernet- 
ted Chain Theory Description of the Dynamics of Ions Con¬ 
fined in Charged Micropores, Phys. Rev. E 59 (1999) 841-849. 
doi:10.1103/PhysRevE.59.841 


[36] W. Olivares-Rivas, P. J. Colmenares, Scaling of Langevin 

and Molecular Dynamics Persistence Times of Inho¬ 
mogeneous Fluids, Phys. Rev. E 85 (2012) 0111171. 

doi:10.1103/PhysRevE.85.011117 

[37] S. Redner, A Guide to First-Passage Problems., Cambridge Uni¬ 
versity Press, Cambridge, United Kingdom, 2001. 

[38] B. Dybiec, E. Gudowska-Nowac, P. Hanggi, Levy- 

Brownian Motion on Finite Intervals: Mean First Pas¬ 
sage Time Analysis, Phys. Rev. E 73 (2006) 046104. 

doi:10.1103/PhysRevE.73.046104 

[39] P. Hanggi, P. Talkner, M. Borkovec, Reaction Rate Theory: 
Fifty Years After Kramers, Rev. Mod. Phys. 62 (1990) 251- 
341. doi:10.1103/RevModPhys.62.251 

[40] B. Jonsson, H. Wennestrom, Computer simulation studies of the 
electrical double layer, in: S. H. Chen, R. Rajagopalan (Eds.), 
Micellar Solutions and Microemulsions. Structure, Dynamics 
and Statistical Thermodynamics, Springer-Verlag, New York, 
1990, Ch. 1, pp. 51-67. doi:10.1007/978-l-4613-8938-5\_3 

[41] K. Sekimoto, Stochastic Energetics, Lecture Notes in Physics, 
Springer, Heidelberg - Germany, 2010. 

[42] T. Sagawa, M. Ueda, Information Thermodynamics: Maxwell’s 
Demon in Nonequilibrium Dynamics, in: R. Klages, W. Just, 
C. Jarzynski (Eds.), Nonequilibrium Statistical Physics of Small 
Systems. Fluctuations, Relations and Beyond, Review of Non¬ 
linear Dynamics and Complexity, Wiley-VCH Verlag, Weinheim 
- Germany, 2013, Ch. 6, pp. 181-211. 

[43] B. Balakrishnan, Fluctuation-Dissipation Theorems from the 
Generalised Langevin Equation, Pramana 12 (4) (1979) 301- 
315. doi:10.1007/BF02894699 

[44] G. Cicotti, J. P. Ryckaert, On the Derivation of the Generalized 
Langevin Equation for Interacting Brownian Particles, J. Stat. 
Phys. 26 (1) (1981) 73-82. doi:10.1007/BF01106787 


15 


